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CHAPTER 1 
INTRODUCTION 


This report has been written to provide documentation to 
the computer code which numerically implements the theory found in 
[1] for the electromagnetic excitation of a thin, straight wire in 
a rectangular cavity which is excited through a small aperture in 
the cavity wall. Because many of the equations and discussions in 
{1] are crucial to an understanding of this code, it is assumed that 
the user has that document available for ready reference. 

This code has been prepared for use on a CDC 6400 computer 
using the RUN fortran compiler. It is also compatable with CDC's 
FIN compiler. However, for later versions of this compiler it may 
be necessary to suppress the non-fatal error message "ARGUMENT TOO 
, SMALL" which is detected in EXP. This error will naturally occur 
during normal program execution and has no effect upon the program 
output. This error can be suppressed by an appropriate call to 
SYSTEMC in program START. 

In addition, every effort has been made to write the code 
in standard fortran. However, it is expected that a few minor modi- 
fications might be needed for its use on non-CDC machines and with 
other fortran compilers. 

Chapter 2 of this document provides a list and description 


” 


of all input parameters for this code. Similarly, chapter 3 outlines 


1 


the various information which the code returns as output. These two 


chapters should provide the user with the needed information to operate 
the code. In addition, chapter 4 discusses the function of each sub- 
program in the code, and a complete listing of all programs can be 


found in the appendix. 


CHAPTER 2 


INPUT PARAMETERS 


File INPUT 

Note that all of the following parameters are read from 
the file INPUT in the main program START. All input parameters that 
represent lengths (A, B, C, XC, YC, R, ZL, ZU, LBDA, XPP, YPP, RA, 
DELFIX) may be input in any system ‘of units desired, so long as the 
same unit is used throughout. Also note that if LBDA is set to unity, 


it has the effect of normalizing all other lengths in units of wave- 


length. 
N -- Number of pulses to be used in wire current 
expansion 
IGRD -- Integer which indicates connection of wire 


to one or both cavity walls. 
If IGRD equals: 


0, wire is unattached 
-1, wire is attached at z=0 

1, wire is attached at z=c 

2, wire is attached at both z=0 and z=c 


A,B,C -- The dimensions of the cavity in the x, y, z 
directions, respectively (note that cavity must 
be oriented such that A<B) 


XC,YC -- The x-y coordinates of the center of the wire 


R -~ The wire radius 
ZL,ZU  -- The z coordinates of the wire endpoints, with 
aL < 2 


LBDA -- The wavelength (X) associated with the frequency 
of operation 


NSM -~ Maximum number of terms used for the sum S3 
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MAX  ~- Maximum value of m permitted in the sum of 
SUBROUTINE BADSUM (See [1], p. 57-58) 


MAXX ~-- Maximum value of m permitted in the sums of 
subroutines GMl, GM2, and GM3 (See [1], p. 57-58) 


CS,LS -- Convergence criteria for all sums (See [1], p. 57-59) 


ADYAN -- Logical variable which controls the method in which 
the components of G, are computed. If ADYAD equals: 


(.T.) first-order difference techniques are used 
to calculate G, from Gy using equation (2.2) from [1]. 


(.F.}) components of Ge are computed directly 


NWALL -- Indicates the cavity wall perforated by the aperture. 
If NWALL equals: 


1, aperture is in x=0 wall (f=) 
2, aperture is in y=0 wall (fi=?) 
3, aperture is in z=0 wall (f=2) 


XPP,YPP -- The x, and x, coordinates of the centerpoint of 
of the Aperturé, such that a, x a, = ff (see [1], p. 23) 
RA -- An array containing the semi-axes of the elliptical 


aperture in the x) and x, dimensions, respectively 


EMAG, PHASAP, THE, PHI,ANG -- Parameters controlling the incident 
plane wave excitation, where 
Penne 


where a, and k are defined in Figure 1 with 
EMAG = |E\|, THE = 6°, PHI = 9", ANG = Wy. 


If PHASAP = (.F.), phase of incident field is referenced 
to the coordinate origin (R = r) and if PHASAP = (.T.), 
the phase is referenced to the centerpoint of the 


aperture (R=r - r,)- Note that all angles are input 
in degrees. 


FIX -- Logical variable which controls whether or not the 
effects of the cavity walls upon the aperture dipole 
Moments are included (See [1], p. 54). If FIX = (.T.), 
such effects are included. 


Figure 1. 


(a) N 


> 


6¢9*, 9°) 


89+, 44) 
(b) 


Input Parameters for Incident Plane Wave. 


DELFIX -- Distance from aperture at which sums are computed 
in order to interpolate the deleted Green's functions 
; (See [1], p. 68) 


CHAPTER 3 
OUTPUT 


File OUTPUT 

The program automatically will create the print file OUTPUT. 
During normal execution, the input information of the previous chapter 
will be printed, as well as the aperture polarizabilities ao, ot, and 
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a (see [1], p. 17-19). The solution for the wire current and the 
22 


total aperture fields (Ep He H ) are given. For comparison, 


v8 1’ 2 sc. sce 
the exterior short circuit fields (E, ; Hy 


2 a are also given 
(see [1], pp. 20-26). 
In addition, several messages of an informative or warning 


Nature may be printed during execution. These messages which can be 


divided into three basic categories, are described below: 


Informative messages -- 
1. NOTE--W(I) INTERPOLATED THRU I = n 

W(I) is an array containing P(a) for a = (I-1)A (see [1], 
p. 63). For small a, P(a) is interpolated as discussed on p.65 of 
{1]. This message indicates that this approximation has been made 
for the first n values of P(a). To reduce n, increase the size of 


the input parameter MAX. 


Non-Fatal Errors -- 
1. WARNING -- DIMENSION SIZE EXCEEDED -- namel TRUNCATED TO n. TO 


FIX, INCREASE SIZE OF name2 ARRAY(S). 


> 
/ 


This message indicates that the arrays used for temporary 
storage in BADSUM, GMl, GM2 and GM3 are not large enough to sum the 
series to the maximum number of terms specified. The particular 
input parameter (MAX or MAXX, specified by namel) is truncated and 
execution is continued. Note that this truncation could subsequently 
cause either fatal error 5 or 6 listed below. To rectify this problem, 
increase the size of the array(s) specified by name2, such that 
MS1 > max(MAX, MAXX) and MS2 > max (c MAX, c MAXX) where c, = B/A and 
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c, = max(A,B,C)/min(A,B,C). 
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Fatal Errors -- 
1. WARNING--MATRIX SIZE EXCEEDS STORAGE ALLOCATION 
The input parameter N is too large. Decrease N and/or 


increase MS such that MS >N+ 3. 


2. WARNING--REORIENT CAVITY SUCH THAT A.LE.B 
A must be less than or equal to B. Redefine coordinate 


system such that this is true. 


5. WARNING--ERROR IN WIRE END POINTS 


Note that ZL < ZU must be true 


4. WARNING--ALL OR PART OF THE WIRE IS OUTSIDE THE CAVITY, OR UN- 
ATTACHED WIRE HAS END POINT ON CAVITY WALL 
Check to see that the input parameters do not violate any 


of the above conditions. 


~ sotiedintts octane oe 


ee 


S. NOT ENOUGH POINTS FOR INTERPOLATION 


Interpolation of P (a) requires its numerical evaluation at 
at least two points. This is not possible with the present input 


value of MAX. To correct this problem, increase MAX. 


6. WARNING--GMI SUM NOT CONVERGED c & 

The value of MAXX is not sufficient to obtain the specified 
convergence (CS,LS) in GMl, GM2, or GM3. Note that c is the conver- 
gence ratio at the last term and 2 is the number of consecutive times, 
if any, the specified convergence ratio CS has been met. To rectify 
this problem, increase MAXX and/or reduce the severtey of the con- 


vergence specified by CS and LS. 


7. WARNING--SUM 3 NOT CONVERGED c 2% 

The value of NSM is not sufficient to obtain the specified 
convergence (CS,LS) in S3. Note that c and & are defined as in 6 
above. To recify, increase NSM and/or reduce the severity of the 


convergence specified by CS and LS. 


8. WARNING--CAVITY DIMENSION TOO SMALL FOR DPLFIX 

This error is caused by insufficient room inside the cavity 
to perform the interpolation necessary for the evaluation of the 
deleted Green's functions in SUBROUTINE DPLFIX. To correct, either 


decrease the size of DELFIX or let FIX = (.F.)., 


File PUNCH 
The program will in certain situations create the punch file 


PUNCH. This file contains the information from the file INPUT, as 
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well as the solution for the wire current and the fields (E »H oH ). 
Ty 1 2 


The file is formatted according to the PUNCH statements found at the 
end of PROGRAM START. The PUNCH file will be created only if execution 
has not been abnormally halted and if sense switch 1 is "fon". This 


can be accomplished on CDC machines with the SCOPE control card ONSW(1). 
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CHAPTER 4 
DESCRIPTION OF PROGRAM AND SUBPROGRAMS 


The programs and subprograms listed in this chapter are those 
used in the code described in this document. Uniformly throughout the 
following description, it is to be understood that variables (X, Y, Z) 
represent Yor (x,y,z) and that (XP,YP,ZP) represent r' or (x',y',z'). 
Similarly, other program variables (denoted in capital letters) ob- 
viously represent physical variables of the problem (for example, A, B 
and C represent the cavity dimensions a, b, and c). Only when such 
a connection is not obvious will it be explicitly*noted. A complete 


program listing is given in the appendix. 


PROGRAM START 

This program's role is primarily that of input/outpyt operations. 
In addition, it also provides for the checking of consistancy among 
the various input parameters. Finally, it directs program execution 
by calling the basic matrix filling and inverting subprograms. It is 
important to note that if the size of the arrays in blank COMMON are 
changed, the appropriate values of MS, MS1, and MS2 should be changed 
accordingly to insure proper error checking. Sufficient dimension size 


can be obtained by adhering to the algorithms listed in chapter 3 under 


non-fatal error 1 and fatal error l. 


SUBROUTINE FILL 


This subroutine fills the matrix eo (see [1], p.52,63) for 
the case of an unattached wire. The necessary values of P(f8) are 


provided by a call to SUBROUTINE GETUM. 


SUBROUTINES FILL1 and FILL2 

These subroutines fill the matrix Oy for the cases of the wire 
attached to the cavity walls at one or both ends respectively. They 
incorporate the necessary modifications to account for non-zero half- 
pulse and the half-testing functions which are needed at the attached 
ends of the wire (see [1], p.50). As in SUBROUTINE FILL, the necessary 


values of P(8) ([{1], p.63) are provided by calling SUBROUTINE GETUM. - 


SUBROUTINE GETUM 

This subroutine fills the array W(I) with the necessary values 
of P(B) to fill the matrix. It will automatically interpolate P(f) 
for sufficiently small 8 (see [1], p.65) using a three-point fit to a 


second order polynominal in 8. 


SUBROUTINE MORFIL 

This subroutine fills the matrices e and er according to (4.7) 
and (4.8) in [1]. In addition it also initializes to i, the identity 
dyad. The necessary components of the dyads have been calculated in 
SUBROUTINE EAP (through calls to SUBROUTINES EDIPL and MDIPL) and are 
contained in the arrays FP, Gl and G2. Note that advantage is taken 
of the fact that G(r") = ae sr) and 8, (F.t") = Mes st) {see le 


Table 2 and (2.7)}. 
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SUBROUTINE INVERS 
This subroutine numerically solves the resulting matrix 


equation (4.5) in [1] using partial-pivoting Gauss Elimination. 


SUBROUTINE FIXUP 

This subroutine fills the matrix o according to (4.9) in [1]. 
This routine is called only if the input parameter FIX = (.T.). The 
values of the deleted Green's functions, which are calculated by 
SUBROUTINE DPLFIX, are supplied to this subprogram through the common 


block/WHYNOT/. 


SUBROUTINE DPLFIX 

This subroutine calculates the deleted Green's functions as 
outlined on pp. 67-70 of [1]. It is called by START only if FIX = (.T.). 
Figures 2 (a and b) show arrangement of points at which the components 
of G, and By are calculated, respectively. Note that A in Figure 2 
represents the input parameter DELFIX. These computations are somewhat 
complicated by the fact that if the aperture is sufficiently close to 
one or more of the side walls of the cavity the choice of points in 
Figure 2 must be modified such that all points remain on the cavity 
interior surface. In addition, when the aperture is close to a side 
wall, it is necessary to remove the image source in that wall before 
interpolation, and then to subsequently add it to the interpolated re- 
sult. Because of these two complications, extensive logic is necessary 


in DPLFIX to adequately handle them. 
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SUBROUTINE EAP 

This subroutine calculates several quantities relating to the 
aperture. It calls subroutines to compute the aperture polarizabilities, 
and the components of G, and 8 necessary to fill O° and o. In addition 
it calls routines which convert the exterior incident field from the 
input parameters (shown in Figure 1) to the x,y,z components of E and H. 


These values are then used to fill the excitation vector in (4.5) of 


{1] into the array E(I). 


SUBROUTINE INC 
This subroutine converts the input specifications for the 
incident field (EMAG,ANG,THE,PHI) to the corresponding components of 


E and H in the spherical polar coordinate system. 


SUBROUTINE CONV 
This subroutine converts a vector from its polar components 
to its cartesian components. It is assumed that the vector is per- 


pendicular to the polar unit vector ?. 


SUBROUTINE ALPHA 
This subroutine computes the aperture polarizabilities for 
the elliptic aperture using (2.15) of [1]. It calls SUBROUTINE ELLIP1iC 


to provide the necessary values of the elliptic integrals. 


SUBROUTINE ELLIPTIC 

This subroutine computes the complete elliptic integrals of the 
first and second kinds of argument X which are denoted EK and E£, re- 
spectively. They are computed using the polynominal approximations of 


2], pp. 591-592. 
[ pp 15 


SUBROUTINE MDIPL 


This subroutine calculates Gh (x, 5¥,) for use in the matrix 


Jz 
e. It fills the array E with these values for p = 1, N. 


SUBROUTINE EDIPL 

This subroutine calculates Coe Tate for use in matrix ig 
It fills the array E with these values for p = 1, N. Note that if 
ADYAD = (.T.), these values are computed indirectly by first computing 


appropriate values of G and using first-order finite difference 


Anz 


techniques. 


FUNCTION S3 


This function numerically computes the function S., defined 


3? 
by (3.2c) in [1], using (3.16) in [1]. 
SUBROUTINES GA, GF, GM and GE 

These subroutines compute the various components of the dyadic 
Green's functions needed in the matrices =. oO and oO (see [1], pp.52-54). 
Note that these subroutines automatically analytically reduce the sum 
that will result in the most rapidly converging remaining double sum 
(see [1], p. 55). Note that IN, SX and SY are dummy arrays and the 
answer is returned as SD. Table l shows the way to call these subroutines 


to compute any desired dyad component. 


SUBROUTINE BADSUM 
This subroutine computes the indefinite integral of the reduced 


kernel P(8) defined by (4.22) in [1]. The m-n plane is subdivided into 
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Table 1. Appropriate subroutine call and value of GEZ to compute 


the various dyad components. 


To Compute 


CALL 
GA(IN,SX,SY,Y, YP,Z2,ZP,X,XP,B,C,A,SD) 
GACIN,SX,SY,Z,ZP,X,XP,%,YP,C,A,B,SD) 
GA(IN, SX,SY,X,XP,Y,YP,Z,ZP,A,B,C, SD) 
GF(IN,SX,SY,X,XP,Y,YP,Z,ZP,A,B,C,SD) 
GF(IN,SX,SY,Y,YP,Z,ZP,X,XP,B,C,A,SD) 
GF(IN,SX,SY,Z,ZP,X,XP,Y,YP,C,A,3,SD) 
GA(IN,SX,SY,Y,YP,2,ZP,X,XP,B,C,A,SD) 
GE(IN,SX,SY,XP,X,ZP,Z,YP,Y,A,C,B,SD) 
GE(IN,SX,SY,XP,X,YP,Y,ZP,Z,A,B,C,SD) 
GECIN,SX,SY,X,XP,Z,ZP,Y,YP,A,C,B,SD) 
GACIN,SX,SY,Z,ZP,X,XP,Y,YP,C,A,B,SD) 
GE(IN,SX,SY,YP,Y,XP,X,ZP,Z,B,A,C,SD) 
GE(IN,SX,SY,X,XP,Y,YP,Z,ZP,A,B,C,SD) 
GE(IN,SX,SY,Y,YP,X,XP,Z,ZP,B,A4,C,SD) 
GA(IN,SX,SY,X,XP,¥,YP,Z,ZP,A,B,C,SD) 
GM(IN,SX,SY,XP,X,ZP,Z,YP,Y,A,C,B,SD) 
GM(IN,SX,SY,XP,X,YP,Y,ZP,Z,A,B,C,SD) 
GM(IN,SX,SY,YP,Y,ZP,Z2,XP,X,B,C,A,SD) 
GUCIN, SX,SY,YP,Y,XP,X,ZP,Z,B,A,C,SD) 
GM(IN,SX,SY,Z?,2,YP,Y,XP,X,C,3,A,SD) 


GM(IN,SX,SY,ZP,Z,XP,X,YP,Y,C,A,3B,SD) 


1T 


GEZ 
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Table 1.~-Continued 


To Compute CALL GEZ 
-8, (r,r') GM(IN,SX,SY,Y,YP,Z,ZP,X,XP,B,C,A,SD) F 
29 Meee! 
g, (@3") GM(IN,SX,SY,Z,ZP,Y,YP,X,XP,C,A,B,SD) F 
XZ 
g, GF") GM, IN,SX,SY,X,XP,Z,ZP,Y,YP,A,C,B,SD) F 
a 
“9, (iF) GM(IN,SX,SY,Z,ZP,X,XP,Y,YP,C,A,B,SD) F 
Vo" 
2, (.F") GM(IN,SX,SY,X,XP,Y,YP,Z,ZP,A,B,C,SD) F 
zx 
ge (r,r') GM(IN,SX,SY,Y,YP,X,XP, Z, ZP,B,A,C,SD) F 
zy 
18 


elliptical bands as discussed on pp. 55-59 of [1] and the double sum ‘ 
is evaluated using the convergence criteria discussed on those pages. 

If the sum does not converge in the allotted number of terms, the 

logical variable CIND is set to (.T.) and control is returned to the 


calling subprogram. SX and SY are input arrays defined as 


SX(I) = sink x sink x', k =21 
x x a 
and 
T 
SY(I) = sink ink y' , k ==T1 
(1) ay ne 
The input X represents 8 and P(8) is returned as SD. 
SUBROUTINES GM1, GM2_ and GM3 


These subroutines numerically sum two dimensional series of 


the following form: 


“oe nm non 
GMI = ab u ¥ sinhy Cc 
mn 
n=n 
) 
where 
sinhy 2, sinhy (c-2,) , I=l 
= i ' = 
Win y,sinhy.2 coshy,z » I=2 
Pycoshy 2 < coshy, (c-z,),1=3 
and 


1 , GEZ = (.F.) 


. , GEZ = (.T.) 


nt 2 2 2 2 2 2 
‘ ky oa ke ky + kK, and 3 = ke - k’. Note that the 


for k_ = = 
x «@ 


input arrays SX and SY are defined as 


SX(I) = u and SY(I) = Y (1-148) 


(I-1+n » 


and a, and n, can each be either 0 or 1. This sum is performed in the 
same fashion as SUBROUTINE BADSUM and the same convergence criteria is 
used. Note that if convergence is not achieved, execution is halted 


and a warning message is printed. 


SUBROUTINES NSS, NCC and NSC 


These subroutines fill the array D in the following ways: 


sink x sink. x' , (NSS) 
x x 
= ' 
D(I) r, cosk y cosk y » (NCC) 
+ ! 
k sink, z coskz » (NSC) 


where 


k = = 2 (k= us (I-1) and k_ #2 21. This subroutine is used to 
x: 38 y b 2 € 
appropriately fill the SX and SY arrays for use in BADSUM, GM1, GM2 


or GM3. 


SUBROUTINE NDE 
This subroutine is a dummy routine to call either NSS, NCC 


or NSC, depending upon the value of the integer L. 
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FUNCTIONS COSH and SINH 


These functions compute the hyperbolic cosine and sine, 


respectively, of argument X. 


FUNCTIONS RSS, RCC, RCS and RSINH 


These functions compute the ratios of various combinations of 


hyperbolic functions. They are defined by 


sinh X sinh Y 


RSG,2) 3 ee 


cosh X cosh Y 


RCC(X,Y,Z) = sinh Z 


and 


sinh X 
sinh Y 


. 


RSINH(X,Y) = 
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